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ABSTRACT 

The electrostatic potential that keeps approximate charge neutrality in neu- 
tron star matter is self-consistently introduced into the formalism for rotochemical 
heating presented in a previous paper by Fernandez and Reisenegger. Although 
the new formalism is more rigorous, we show that its observable consequences 
are indistinguishable from those of the previous one, leaving the conclusions of 
the previous paper unchanged. 

Subject headings: stars: neutron — dense matter — relativity — stars: rotation 
— pulsars: general 

1. INTRODUCTION 

Rotochemical heating is one of several mechanisms that can keep neutron stars (NSs) 
hot beyond their standard cooling time of ~ 10 7 yr (see Fernandez & Reisenegger 2005 
[hereafter FR05] for further references, and Yakovlev & Pethick 2004 for a recent review 
on neutron star cooling). It originates in deviations from beta equilibrium produced as the 
star gets compressed due to its decreasing rotation rate, implying a conversion of rotational 
energy into thermal energy. This effect is capable of keeping old millisecond pulsars (MSPs) 
at surface temperatures ~ 10 5 K, consistent with likely thermal ultraviolet emission from 
PSR J0437-4720 (Kargaltsev, Pavlov, & Romani 2004; FR05). 

The heating mechanism works as follows. As a NS spins down, the centrifugal force 
acting on each fluid element diminishes, raising the local pressure. This compression results 
in a change of the equilibrium concentration of each particle species. Reactions that modify 
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the chemical composition drive the system towards the new equilibrium configuration. But if 
the rate at which reactions do this task is slower than the change of the equilibrium concen- 
trations due to spin-down compression, the matter remains out of chemical equilibrium. The 
excess energy stored in the chemical imbalance is dissipated, partly by enhanced neutrino 
emission and partly by heating the matter in the star. Thus, the star is kept warm and can 
radiate for long after the standard cooling time. 

After its introduction (Reisenegger 1995), this heating mechanism was studied by several 
authors: Cheng & Dai (1996) applied it to quark stars, Reisenegger (1997) estimated the 
effects of superfluidity, and Iida & Sato (1997) studied the heating due to compositional 
transitions in the crust due to spin-down compression. 

The most striking prediction associated with rotochemical heating is that, if the spin- 
down timescale is substantially longer than any other timescale involved (with the likely 
exception of magnetic field decay), the star arrives at a quasi-steady state, in which the 
temperature depends only on the current, slowly changing value of VLVL (the product of the 
angular velocity and its time derivative), and not on the star's previous history (Reisenegger 
1995). This provides a simple way to constrain the physics involved in theoretical models, 
once the spin parameters and observed surface temperature of a given, sufficiently old NS 
are known. 

In spite of the previous work, FR05 were the first to go beyond order-of-magnitude 
estimates and calculate the thermal evolution of NSs with rotochemical heating, taking the 
structure of the star fully into account in the frame of general relativity, and using realistic 
equations of state of dense matter. They developed a general formalism to treat the evolution 
of the temperature and departures of chemical equilibrium, implementing it for the simplest 
possible core composition, namely neutrons, protons, electrons, and muons (npefi matter). 
Using state-of-the-art equations of state, and modified Urea reactions as the dissipative 
process, they obtained results consistent with the observed ultraviolet emission from the 
millisecond pulsar PSR J0437-4720 (Kargaltsev et al. 2004). 

Here, we show that their formalism is not completely rigorous, in the sense of ignoring 
the perturbations to the electrostatic potential that must be present inside the NS in order 
to keep nearly exact charge neutrality among the different charged particles present inside 
the star. 

We start (in §2) by discussing the physical origin, importance, and strength of the 
electrostatic potential expected inside a NS, as well as its relation to the (nearly exact) 
condition of local charge neutrality. Then (in §3), we review the main components of the 
rotochemical heating formalism developed by FR05, which are rigorously correct and not 
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directly affected by the presence of electrostatic potential perturbations, showing that their 
main result, namely the existence and properties of the quasi-steady state, are unchanged. 
In §4, we introduce the electrostatic potential into the equations of FR05 that relate the 
different chemical imbalance variables (chemical potential imbalances and deviations from 
equilibrium particle abundances) to each other. The quantitative consequences of these 
modifications are explored in §5, whereas §6 contains our main conclusions. 

2. ELECTROSTATIC POTENTIAL IN A NEUTRON STAR 

It is well-known that a neutron (n) in vacuum is an unstable particle, decaying into a 
proton (p), an electron (e), and an electron antineutrino (P e ). In NSs, neutrons are stabi- 
lized by the presence of enough protons and electrons to block (through the Pauli exclusion 
principle) the possible final states of the decay. Neutron and protons have high masses and 
are therefore held inside the star by its gravitational potential well, whose depth must be 
larger than the kinetic part of their Fermi energies. However, this potential is not deep 
enough to hold the much lighter electrons, whose Fermi energy is relativistic (much larger 
than their mass), and which would therefore tend to escape, yielding a net positive charge 
on the NS. However, at NS densities, already a small charge imbalance between protons and 
electrons causes an electrostatic field that can hold in the electrons and at the same time 
partially balance the gravitational force on the protons, preventing them from sinking into 
the center of the star, due to their relatively small kinetic energies. Thus, a non-rotating NS 
in hydrostatic, chemical (beta), and diffusive equilibrium contains a spherically symmetric 
electrostatic potential ijj(r) that keeps the charged particles near local charge neutrality. The 
existence of such an electric potential in stars has been known to astrophysicists since the 
early twentieth century. A recent discussion, including early references, is found in Neslusan 



More formally, and using a relativistic description of gravity, in contrast to the Newto- 
nian description in the previous paragraph, we describe the star by the spherically symmetric 
metric 



(Schutz 1985). The condition of diffusive equilibrium for each particle species i is that its 
total, redshifted chemical potential, 
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is uniform throughout the star. Here, qi is the charge of each particle species, and /ij is its 
intrinsic chemical potential (Kittel & Kroemer 1980), determined by the kinetic energy of 
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the species (or equivalently, by their number density rii) and their local interactions with 
particles of the same and other species (which are determined by all the densities). For a 
given metric and electrostatic potential, these conditions determine the distribution of the 
individual particle species by fixing the spatial dependence of the intrinsic chemical potential. 
For a typical density inside NSs, n e ~ 10 37 cm" 3 , the extremely relativistic electron chemical 
potential /i e ~ ^c(37r 2 n e ) 1 / 3 10 2 MeV, which sets the scale for the electrostatic potential, 
~ n e /e ~ 10 8 V, where e = q p = —q e . 

The field equation for the electrostatic potential is Maxwell's equation, F .* = 47r J", 
where the components of the electromagnetic field tensor can be written in terms of the 
4-vector potential components A 13 as F af3 = Ap. a — A a .p, and J a are the components of 
the charge-current 4-vector. In the static system of interest here (ignoring currents and 
magnetic fields inside the neutron star), the space components of the 4- vectors can be set to 
zero: A 1 = = J 1 , and the time component A = tp. A substantial amount of algebra yields 
the relativistic Poisson equation for the electrostatic potential, 



where R a p are the components of the Riemann curvature tensor, and the primes indicate 
derivatives with respect to the radial coordinate r. The terms on the left-hand side (the 
generalization of the laplacian operator of flat spacetime) are of order \ip\/R 2 , where R ~ 
10 6 cm is the radius of the NS. Thus, a minute excess positive charge density ~ fi e /(eR) 2 ~ 
10 2 cm" 3 (compared to the typical density of each charged species, ~ 10 37 cm 3 mentioned 
above) is enough to produce the required potential. Therefore, to an excellent approximation, 
we can assume that the electrostatic potential is exactly that required to keep local neutrality 
everywhere among the charged particle species, replacing equation (3) by the condition of 
charge neutrality 



imposed on the densities and obtaining the non-zero electrostatic potential from this 
condition, together with equation (2). 

3. CHEMICAL IMBALANCE AND ROTOCHEMICAL HEATING 

The formalism of FR05 can be summarized as follows (we refer to the original paper 
for more details). The changing rotation rate Q(t) of the star, assumed to be much slower 
than the "break-up" or "Keplerian" rate, slightly perturbs its structure, making the total, 
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chemical-equilibrium number N^ q of particles of species % in the star change with time as 

N e t q = 2I nt Sl(l, (5) 

where Iq^ is a constant, depending on the structure of the unperturbed, non- rotating star. 
The actual number of particles of the same species, iVj, is changed by non-equilibrium weak 
interactions in the stellar core as 

N i= f dVe*J2 Ar *- ( 6 ) 

J core a 

The integral is performed over the volume of the unperturbed, non-rotating stellar core, and 
the sum runs over different kinds of reactions (labeled by a), each of which creates a net 
number AI^ of particles of species % per unit time and per unit volume. The latter quantity 
is a function of the local density, temperature, and, most importantly, of the local chemical 
imbalance, quantified by the variables i] np i = fi n — fi p — fii, where I stands for the charged 
leptons: electrons (I = e) or muons (/ = //). If the r] np i are non-zero, so are the net reaction 
rates, with their signs such that the chemical imbalances tend to be decreased. The evolution 
of the excess (or deficit) of particles of species i, 5Ni = N { — iVf 3 , is therefore determined by 
the competition between the changing rotation rate of the star, which tends to increase its 
absolute value, and the reactions, tending to decrease it: 




The evolutionary time scale of the star is assumed to be much longer than the time scales 
for diffusion of particles and heat across it. So, in particular, the redshifted temperature 
rpoo _ y( r ) e *(»0 i s considered to be uniform inside the star, except for a thin surface layer. 
This temperature evolves as 

T°° = i [US - L~ - L~] , (8) 

where C is the total heat capacity of the star, is the total power released by reactions 
and other heating mechanisms, is the total neutrino luminosity, and is the photon 
luminosity, which are functions of the chemical imbalances i] np i and the temperature T 
everywhere inside the star. 

Although the particles are in diffusive equilibrium, moving easily everywhere inside 
the star, they are generally not in chemical equilibrium (free to convert into one another). 
Therefore, non-zero, but spatially uniform, redshifted chemical imbalances are set up, 

Vn P i =rf?-v£-ri° = Vn P ie", (9) 
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not involving the electrostatic potential, which cancels out due to the equal and opposite 
charges of protons and leptons. Being uniform in the star, these redshifted chemical im- 
balances are the ideal variables to quantify the departure from chemical equilibrium and 
follow its time evolution. The reaction rates on the right-hand side of equations (6) and (7) 
can be written in terms of rf^ pl and T°° instead of the corresponding un-redshifted, local, 
position-dependent variables. 

In order to close the system of dynamical equations, it is necessary to write the chemical 
imbalances, 77^, in terms of the excess numbers of particles, SNi, which is where the inac- 
curacy of the previous treatment (FR05) is located, as we shall discuss below. However, at 
this point, we may note that, if we are only interested in the quasi-steady state reached by 
NSs with slow spin-down evolution, we only need to set the right-hand sides of equations (7) 
and (8) equal to zero, immediately yielding a non-degenerate set of algebraic equations for 
r}™ t and T°°, which is equivalent to that of FR05 and therefore has identically the same 
solutions. 



4. RELATIONS AMONG THE VARIABLES QUANTIFYING THE 

CHEMICAL IMBALANCE 

When the chemical equilibrium is slightly perturbed, the total, redshifted chemical 
potential of equation (2) is displaced from its equilibrium value by the small amount 

5nT = [SfU + qM + (tu + g^)5$]e*, (10) 

where the prefix 5 refers to the respective variable being slightly perturbed with respect to 
the chemical equilibrium state, maintaining diffusive equilibrium (so Sfj,f° remains uniform). 
We note that equation (9) of FR05 considered only the first term. If some neutrons are 
converted into proton-electron pairs (plus escaping antineutrinos), its effect on the energy 
density and pressure, and therefore on the structure of the star and the spacetime metric, 
is not large, as verified explicitly (in the Newtonian approximation) in Appendix A. We 
therefore invoke the relativistic Cowling approximation (e. g., Finn 1988), setting 5<E> = 0. 
However, the absolute value of the electrostatic potential must increase, in order to hold 
more protons and electrons together than before. This effect is considered explicitly in what 
follows. 

The intrinsic chemical potential perturbations are related to local density changes, Srii = 
Y2i(dni/d[/,j)6[/,j, where, from equation (10), SfXj = <5/i°°e~ <I> — qjSip. Imposing conservation 
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of local charge neutrality, J2i Qi^ n i — 0, yields 



<^ = ^^ 775 /a \ • 11 

Lr, S 9rg S (<9n r /(9/X s ) 

Thus, after some algebra, the perturbation to the total number of particles of one species in 
the star can be written as in equation (11) of FR05: 

SNi = J dVSm = Y^BijSfi?, (12) 



but with the coefficients redefined as 

drii Y.km ( lkqm{.dn i /dix m ){dn k /dix j ) 



= I dVe 
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(The second term inside the square brackets is missing in eq. [12] of FR05.) The thermo- 
dynamic identity drii/dfij = drij/d/ii can be used to show that B^ = Bji. More impor- 
tant, however, is to note that the 4 x 4-matrix B is singular, since charge conservation is 
"built in," in the sense that, for any combination of four values for the variables SfJ-f, it 
is guaranteed that 'Yl li qi5N i = 0, i. e., three of the rows of B are not linearly independent: 
B p j = B e j + B^j. One manifestation of this singularity is the existence of "trivial" per- 
turbations with <5/i°° = qiSipe® , for which 5 = and therefore the particle distributions 
inside the star are not changed. These perturbations can be interpreted as redistributions 
of charges outside the star, or just as redefinitions of the zero point of the potentials. 

In order to obtain an invertible problem, we reduce our four degrees of freedom to 
three, by defining the new variables 5jl™ = <5/i£° and 5fi^° = S/i™ + 5fj,f°, the latter for 
both / = e, //, and eliminating the row and column corresponding to the protons from 
the matrix B, in this way converting it into a symmetric and invertible 3 x 3-matrix B, 
with 5Ni = J2j= n e a Bij^f for i — n,e, /i. (The proton density and total number can be 
recovered from those of electrons and muons through the condition of charge neutrality.) 
The components of the inverse matrix, B _1 , relate the chemical imbalance variables to the 
particle number perturbations: r)™ pl = 6p% - 5fi™ = E i=n , e# (^ 1 - B u l )5Ni. 

Finally, this can be further reduced by considering that the reactions also conserve 
baryon number 1 , i. e., SN n = —8N P = —5N e — 5N^, which allows to also eliminate SN n: 
writing 

Vn pl = E " A7 1 - Kn + B£)6Ni. (14) 

i=e,fi 



1 Note that, contrary to charge neutrality, which can be imposed locally everywhere in the star, baryon 
number is only a globally conserved quantity. 
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We can make contact with the notation of FR05 by writing 

V ™ e = -Z npe 5N e - Z np SN^ and rj^ = -Z np 5N e - Z npfl SN^ (15) 

with the coefficients rewritten as 

Z n P = B^-B^-B^ + B^, (16) 

Z npe = B~l - 2B~l + B'e, (17) 

Z npf i = B n l — 2B n ^ + B^ , (18) 

instead of equations (54) through (56) of FR05. These re-definitions were chosen so that 
their equations (52) and (53) for the time evolution of the chemical imbalance variables rf^L 
remain formally correct. The coefficients W np i, which determine the rate by which chemical 
imbalances are built up by the changing rotation rate, also remain correctly expressed by 
equations (57) and (58), although more easily written as 

W npe = Z npe l^ e + Z np In^ and W np ^ = Z np lQ e + Z npil I^^. (19) 



5. QUANTITATIVE EVALUATION AND DISCUSSION 

Formally, the only changes that have to be introduced to the equations of FR05 are to 
replace their equations (9), (12), (13), (54), (55), and (56) by the corresponding equations 
in the previous sections. In their dynamical system for the time evolution of the variables 
T°° and 77^, the only difference lies in the values of the coefficients Z np , Z npe) and Z npfl , 
now defined by our equations (16), (17), and (18), instead of their equations (54), (55), and 
(56), and in the coefficients W np i, which depend on the former (e. g., our eq. [19]). We 
have argued in §3 that the quasi-steady state eventually reached by old NSs is unaffected 
by these changes. However, the rate at which this state is approached (directly depending 
on the coefficients W np i\ see, e. g., FR05, eqs. [76] and [81]) could be affected, as well as the 
less intuitive time evolution of younger NSs with stronger magnetic fields (e. g., Reisenegger 
1995). 
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Fig. 1. — Coefficients Z np , Z npe , and ^ npAt as functions of stellar mass, for models built with 
the A18 + 5v + UIX* equation of state (Akmal et al. 1998). Dashed lines correspond to the 
formalism of FR05 (without explicit account of electrostatic potentials), while solid lines are 
the results of the corrected formalism presented in the present paper. 
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Figure 1 shows a comparison of the coefficients Z np , Z npe , and Z nm , as functions of 
stellar mass, for one particular equation of state (Al8+<5-u+UIX* , Akmal et al. 1998), for 
the formalism of FR05 (dashed lines) and the more rigorous one presented here (solid lines). 
It can be seen that Z np hardly changes, whereas Z npe and Z npfJi show changes whose fractional 
amplitude increases with mass, up to a ~ 10 — 20 % level at the largest masses. Calculations 
with other equations of state yielded very similar results. 
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Fig. 2. — Evolution of the internal redshifted temperature, T°°, and chemical imbalance 
variables rf^Jk and rj^/k, for NS models of 1.47 M (left) and 2.13 M© (right), calculated 
with the A18 + 5v + UIX* equation of state (Akmal et al. 1998), with initial conditions of 
vanishing chemical imbalances and = 10 8 K at t — 0, and assuming magnetic dipole spin- 
down with field strength B = 10 8 G and initial period Pq — 1 ms. Dashed lines represent 
results for the formalism of FR05, whereas solid lines correspond to the present, corrected 
formalism. 
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In Figure 2, we compare the time evolution of the chemical imbalance variables, ?7^ e 
and rf£ , as well as the redshifted internal temperature, T°°, as obtained in both formalisms 
for two different neutron star masses. For the lower mass, particles inside the star undergo 
only (slow) modified Urea reactions, whereas at the higher mass electrons can undergo (fast) 
direct Urea processes, whereas muons are still restricted to modified Urea. This accounts 
for the relatively simple evolutionary curves on the left, as compared to the more complex 
one on the right. (See FR05 for a more detailed discussion.) It can be seen that, consistent 
with our previous discussion, there is no evidence for differences between the two formalisms 
in the asymptotic, quasi-steady regime, whereas some differences are seen at the epochs of 
strong time evolution, especially in the higher mass star. 

6. CONCLUSIONS 

We have shown that the electrostatic potential inside a NS plays an important role in 
determining the spatial distribution of charged particles and their nearly exact local charge 
neutrality everywhere in the star. We take it into account explicitly in an extension of the 
formalism for rotochemical heating of FR05, in which the relations between the two sets 
of variables quantifying chemical imbalances (chemical potential differences and departures 
from equilibrium abundances) are affected, without changing the rest of the formalism. The 
quasi-steady state reached by old NSs is strictly unaffected by these corrections, and the 
time-dependent evolution at earlier times shows quantitatively negligible changes. 
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we address in the present paper. We also thank Olivier Espinosa for useful discussions, 
and Juan Veliz for assistance with the preparation of figures. This work received financial 
support from FONDECYT Regular Grants 1020840 and 1060644, from the FONDAP Center 
for Astrophysics, from the Visiting Professors Program of Pontificia Universidad Catolica de 
Chile, from the Russian Leading Scientific School grant 9879.2006.2, and from the Russian 
Fund for Basic Research grant 04-02-17590. 



A. Quantitative evaluation of gravitational potential perturbations 



-13- 



1 




i 1 , 1 1 1 , 1 1 1 1 1 , 1 1 1 

2 4 6 8 10 12 14 

R [km] 

Fig. 3. — Perturbation of the gravitational potential in a neutron star model with the first 
PAL equation of state (Prakash et al. 1988) and total mass M — 1.6 M , subject to a 
chemical perturbation corresponding to a 1% increase of the electron and proton density 
everywhere in the star, together with a corresponding decrease in the neutron density, so as 
to conserve total baryon number. 
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Here we verify whether the gravitational potential (or metric) perturbation in equation 
(2) can indeed be neglected. We consider a neutron star model of 1.6 M & , calculated in full 
General Relativity with the first PAL equation of state (Prakash et al. 1988). For simplicity 
using perturbation theory in a Newtonian framework, we calculate the gravitational potential 
perturbation <5$ corresponding to an arbitrarily chosen deviation from chemical equilibrium, 
taken as an excess electron and proton density of 1% everywhere (muons are not present in 
this model), with a corresponding decrease in the neutron density, so as to conserve baryon 
number. Figure 3 shows 5$ as a function of the radial coordinate. For comparison, 5fj, e / ' ji e ~ 
5n e /(3n e ) rs 0.003, so the intrinsic chemical potential perturbation term in equation (2) is 
about 50 times larger than that corresponding to the gravitational potential. 
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